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Abstract 

We discuss several phenomenological aspects of the fully constrained version of the 
next-to-minimal supersymmetric extension of the standard model (cNMSSM). Assum- 
ing universal boundary conditions at a high energy scale for the soft supersymmetry- 
breaking gaugino, sfermion and Higgs masses as well as for the trilinear interactions, 
we find that the model can satisfy all present constraints from colliders and cosmo- 
logical data on dark matter, B- and muon-physics. The phenomenologically viable 
region of the parameter space of the cNMSSM can be described by essentially one 
single parameter as the universal gaugino mass parameter Mi/27 corresponds to 
small values for the universal scalar mass niQ. The lightest supersymmetric particle is 
always a singlino-like neutralino that is almost degenerate with the lightest tau slep- 
ton. We study the particle spectrum of the model and its signatures at the LHC, such 
as possibly long-lived tau sleptons at the end of decay chains, that would distinguish 
the cNMSSM from the constrained MSSM. 
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1 Introduction 



The next-to-minimal supersymmetric standard model (NMSSM) [1-3], in which the spec- 
trum of the minimal extension (MSSM) is extended by one singlet superfield, was among 
the first supersymmetric (SUSY) models based on supergravity-induced SUSY-breaking 
terms. It has gained a renewed interest in the last decade, since it solves in a natural and 
elegant way the so-called fi problem [4] of the MSSM: in the NMSSM, this parameter is 
linked to the vacuum expectation value (vev) of the singlet Higgs field, and a value close 
to the SUSY-breaking scale is dynamically generated. In fact, the NMSSM is the simplest 
supersymmetric generalization of the standard model (SM) in which the SUSY breaking 
scale is the only scale in the Lagrangian, since it allows for a scale invariant superpotential. 

In contrast to the non or partially constrained versions of the NMSSM that have been 
intensively studied in recent years [5] and which involve many free parameters, the con- 
strained model (cNMSSM) has soft SUSY-breaking parameters that are universal at a high 
scale such as the grand unification (GUT) scale: common gaugino (M1/2) and scalar (mo) 
masses as well as trilinear couplings (Aq), as motivated by schemes for SUSY-breaking that 
are mediated by flavour blind gravitational interactions. If one assumes that the soft SUSY- 
breaking parameters involving the additional singlet sector are universal as well, the model 
has the same number of unknown parameters as the constrained MSSM (cMSSM) [6]. 

General features of the constrained NMSSM parameter space as well as aspects of its 
phenomenology have been discussed earlier in Refs. [2,3]. These studies already revealed 
that the allowed range for the parameters M1/2, mo and Aq is different from that of the 
cMSSM. Small values for mo are disfavored in the cMSSM, as they lead to charged sleptons 
that are lighter than the neutralino Xi? the preferred lightest SUSY particle (LSP). In the 
cNMSSM, small mo is needed to generate a non- vanishing vev of the singlet Higgs field [2]. 
The slepton LSP problem can be evaded owing to the presence of the additional singlino-like 
neutralino which, in large regions of the cNMSSM parameter space, is the true LSP [7]. 

Since the early studies of the cNMSSM, bounds on the Higgs and SUSY particle spec- 
trum from high-energy collider data and low-energy measurements have become more se- 
vere, and the dark matter relic density has been determined quite accurately [8] . In addition, 
tools for a more precise determination of the mass spectrum and couplings [9], and for the 
computation of the dark matter relic density [10], have become available. In fact, the tools 
which allow to calculate the spectrum in the fully constrained NMSSM have only been de- 
veloped recently and a brief account of the resulting phenomenology has been reported in 
Ref. [11]. In this paper, we investigate in more detail the parameter space of the cNMSSM 
in light of the recent constraints, using the updated tools. 

A priori, it is not obvious if the latest constraints on the Higgs sector from LEP [12] and 
on the dark matter relic density [8] can be simultaneously satisfied in the fully constrained 
cNMSSM with its additional CP-even and CP-odd Higgs states and a singlino-like LSP. By 
scanning the cNMSSM parameter space with with the program NMSSMTools [9], which 
calculates the Higgs and SUSY particle spectra in the NMSSM, we found regions where this 
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is indeed possible [11]. These phenomenologically viable regions correspond to a regime of 
very small or vanishing universal scalar mass mo, and the requirement of a correct relic 
density determines the universal trilinear Aq as a function of the universal gaugino mass 
Ml/2. For low values of M1/2, for which Aq is also quite small, the soft SUSY-breaking 
terms are necessarily close to their "no-scale" values, Ao = mo = [13]. (mo = 0, but 
nonvanishing values for M1/2 and Aq at a high scale, can also originate from a strongly 
interacting conformal hidden sector [14]). For small enough values of the Yukawa coupling 
between the two doublet and the singlet Higgs fields, A < 10~^, constraints from LEP on 
the Higgs sector as well as constraints from i?-physics are also satisfied. Notably, the model 
is very predictive in the sense that the complete sparticle spectrum depends essentially on 
just one parameter, which can be taken as M1/2. 

Moreover, in these viable regions of the parameter space, the cNMSSM can also explain 
the deviation of the observed anomalous magnetic moment of the muon {g — 2)^ with 
respect to the SM expectation [15] as well as the 1.7a and 2.3 a excesses of events observed 
at LEP [12, 16] corresponding to Higgs masses around 115 GeV and 98 GeV, respectively. 

Some of these results have already been published in Ref. [11]; here we discuss in more 
detail the full allowed parameter space, the complete sparticle and Higgs spectrum, as well 
as the phenomenological implications for future experiments, notably for the LHC. 

An unconventional but general property of the sparticle spectrum of the cNMSSM is a 
singlino-like LSP with very small couplings to non-singlet particles, and a stau next-to-LSP 
(NLSP) with a mass close to that of the LSP. The stau will appear in practically all sparticle 
decay cascades. The small value of mo as compared to M1/2 implies that all squarks are 
lighter than the gluino, a feature which will be relevant for sparticle searches at the LHC. 
The SM-like CP-even Higgs boson has a mass in the 115-120 GeV range. 

In some regions of the parameter space of the cNMSSM, for instance for a very small 
Yukawa coupling A, the stau lifetime can be very large, with visibly displaced vertices 
originating from its decay. For larger values of A, the additional mostly singlet-like CP- 
even Higgs boson can have sizable couplings to gauge bosons and a mass around 98 GeV, 
which would in turn explain the 2.3 a excess of events observed at LEP [12, 16]; the other 
1.7 cr excess of Higgs-like events would be due to the nearly SM-like next-to-lightest CP-even 
Higgs with a mass close to 115 GeV. 

Hence, not only the sparticle spectrum of the cNMSSM should already allow to discrim- 
inate it from the cMSSM at the LHC, but additional unconventional phenomena such as 
displaced vertices or a more complicated Higgs sector can also occur. 

The layout of the paper is as follows: in Section 2 we describe the model, its free 
parameters and discuss the phenomenologically viable cNMSSM parameter space with the 
help of analytic approximations. In Section 3, we present the results for the Higgs and 
sparticle spectra. In Section 4 we discuss additional phenomenological aspects of the model 
that are relevant to the LHC: sparticle production cross sections and decay cascades, the 
possibility of displaced vertices from stau decays and features of the Higgs sector. Finally 
we comment on tests at e~^e~ colliders and the direct and indirect detection of dark matter, 
which can rule out the present model. A summary and concluding remarks are presented 
in Section 5. 
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2 The constrained NMSSM 



2.1 The parameters of the cNMSSM 

We consider the NMSSM with a scale invariant superpotential given by 

W = htQHj% - hQH^b^ - KLHar^ + XSHuHa + | , (2-1) 

where hatted letters denote superfields. H^, and S represent the complex scalar Higgs 
fields, with hd and s their vacuum expectation values. Tilded letters will denote the 
scalar components of quark and lepton superfields. For simplicity, only third generation 
(s)fermions have been included and Q,L stand for the (t, 6) and (r, z/,-) SU(2) doublet 
superfields. The three first terms in eq. (12. ip are the usual generalization of the Yukawa 
interactions, while the last two terms involving the singlet superfield S substitute the iiHuHd 
term in the MSSM superpotential: a non-vanishing value s at the minimum of the Higgs 
potential generates an effective /i term 

/^eff^As. (2.2) 

The singlet superfield S contains a neutral CP-even and a neutral CP-odd scalar, as 
well as a neutralino. All these states mix with the corresponding components of the Hu and 
Hd superfields, increasing the rank of the CP-even, CP-odd and neutralino mass matrices 
by one as compared to the MSSM. Conventions for signs and mixing matrices are chosen 
as in the SLHA2 [17] convention, and we take A > 0. 

The soft SUSY-breaking terms consist of mass terms for the gaugino, Higgs and sfermion 
fields (for the latter, we will use the notation of the third generation; a sum over the three 
generations is implicitly assumed) 
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Id 

+ ml\bl\+ml\L'\+ml\fl\ , (2.4) 

as well as trilinear interactions between the sfermion and the Higgs fields, including the 
singlet field 

-Aril = (^htAtQHui%-hbAbQHdb%-hrArLHdf^ 

+ XAxHuHdS + ^kA^ S'^ +h.c. . (2.5) 

All parameters in the above Lagrangian depend on the energy scale via the corresponding 
RG equations, so that the dominant radiative corrections involving large logarithms are 
accounted for. In the fully constrained cNMSSM, one imposes unification of the soft SUSY- 
breaking gaugino masses, sfermion and Higgs masses as well as trilinear couplings at the 
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grand unification scale Mqut^ 

Ml = M2 = M3 = Ml/2 , 

rriH^ = niHa = mg = mq = rut = nib = uil = rUr = mo , 

At = Ab = Ar = Ax = A, = Ao. (2.6) 

Then, apart from gauge and quark/lepton Yukawa couplings, the Lagrangian of the 
cNMSSM depends on the five input parameters 

Ml/2 , mQ , Aq , X and k . (2.7) 

Requiring the correct value of Mz reduces the dimension of the parameter space from five 
to four. 

In principle, one could start with four independent parameters (such as mo/Mi/2, 
Aq/Mi/2, a and n), and integrate the RG equations for all soft terms from Mqut down 
to the SUSY scale Msusy, defined by the order of magnitude of the soft SUSY-breaking 
terms. Subsequently, one could minimize the effective potential with respect to hu, hd and 
s and determine the overall scale of the soft terms in eq. ( 12. 7p from the correct value of Mz, 
as done in Ref. [2]. However, since tan/5 = h^/hd is then obtained as output (while the 
top quark Yukawa coupling ht would be an input), it becomes very difficult to obtain the 
correct value for the top quark mass rrit or, given mf, ht can only be obtained once tan/5 is 
known. Since ht is very important for the radiative corrections and the RG evolution, it is 
much more convenient to allow for tan (3 as an input parameter, which permits to determine 
ht at the weak scale from the beginning in terms of m^. 

All in all, the following procedure is feasible in practice: apart from Mz-, the five pa- 
rameters 

Ml/2 , mo , Ao , X and tan/5 (2.8) 

are allowed as inputs. The parameters k, the soft singlet mass m| as well as the vev |s| 
(or |/ieff| = A|s|) are determined at MgusY through the three minimization equations of 
the scalar potential with respect to hu, hd and s. (With the convention X > 0, k typically 
turns out to be positive as well, and of C(A/10); the sign of s or can still be chosen at 
will.) This is the procedure employed by the routine NMSPEC within NMSSMTools [9], 
which calculates the spectra of the Higgs and SUSY particles in the NMSSM in terms of 
the soft SUSY breaking terms at Mqut (except for the parameter m|), tan/5 at the weak 
scale (defined by Mz) and A at the SUSY scale Msusy- 

Clearly, the soft singlet mass squared m| at Mgut will not coincide with ttiq in general 
(for a recent analysis allowing for a non-universal singlet mass term, see Ref. [18]). However, 
one can confine oneself to regions in parameter space where the difference between m| and 
TTT-Q is negligibly small. This condition leaves us with an effective 4-dimensional parameter 
space, consistent with the considerations above. In practice, we determine tan/5 by the 
requirement that m| at Mqut should be close to ml: we impose |m|(MGUT) — ml\ < 
(5 GeV)^, which typically requires to tune the fourth decimal of tan/5. This should not be 
interpreted as a fine-tuning, since m| should be considered as an input parameter, whereas 
tan (3 is determined by the minimization of the effective potential. 

For the most relevant SM parameters, the strong coupling and the bottom/top quark 
masses, we chose a^iMz) = 0.1172, mp(mb) = 4.214 GeV and mf^^ = 171.4 GeV [19]. 
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2.2 Constraints from the scalar potential 

Let us begin by recalling some conditions on the parameters M1/2, ""^o and Aq of the 
cNMSSM, which follow from a phenomenologically acceptable minimum of the Higgs po- 
tential [2]. First, the vev s of the singlet has to be non- vanishing. The dominant s-dependent 
terms in the Higgs potential are given by 

V{s)r^ K^s^ + lKA^s'^ + mls^ + ... , (2.9) 

and one easily finds that the condition for a non-vanishing value for s at the absolute 
minimum is equivalent to the inequality 

ml < . (2.10) 

For small A and k (as will be the case below), the parameters and rris are hardly 
renormalized between the GUT and the electroweak scales, and the above condition trans- 
lates into the first constraint on the parameter space (assuming rriQ > 0) 

mo<l\Ao\. (2.11) 

Next, we consider the CP-odd Higgs boson mass matrix. The dominant term in its 
singlet-like diagonal element is given by [1-3, 9, 17] 

-SkA^s, (2.12) 

and must be positive. For positive s and k this implies negative trilinear couplings 

A,r^Ao<0. (2.13) 

We must also consider the constraints on the parameter space arising from vacuum 
stability. Dangerous instabilities of the scalar potential along charge and colour breaking 
(CCB) directions in field space can occur [2,20], notably for large values of |y4o|. The most 
dangerous CCB direction is along the D-flat direction 

\E,\ = \L^\ = \H,\ (2.14) 

in field space, where Ei and Li are the right- and left-handed selectron fields, and where 
the term ~ hgAe in the scalar potential can give a large negative contribution. Once Ag 
and the corresponding soft masses at the appropriate scale (using the corresponding RG 
equations) are expressed in terms of Aq, M1/2 and mg, the condition for the absence of such 
a charge and colour breaking minimum becomes [2] 

(A0-O.5M1/2)' < 9m2 + 2.67M2/2 . (2.15) 

In our analysis we will obtain relatively small values for Aq {Aq ~ —^Mi/2, cf. Fig. [1] 
below), for which eq. fl2.15p is satisfied independently of the value of tuq. 
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More delicate could be unbounded-from-below (UFB) directions in field space, which 
are both D-flat and F-flat. In Ref. [20], it has been clarified that such dangerous directions 
in the field space of the MSSM are still present in the NMSSM, although the singlet vev s 
gives an additional positive contribution to the potential. Analytic approximations to the 
potential along such dangerous directions have been studied in Ref. [21], with the following 
results: 

i) the inequality tuq > 0.3Mi/2 (for the large tan/3 which will be relevant here) is an 
approximate condition for the absence of deeper minima in these directions - this 
inequality will be violated below; 

ii) the decay rate of the standard vacuum is usually much larger than the age of the 
universe; hence we have to assume that the early cosmology (temperature-induced 
positive masses squared for the squarks and sleptons) places us into the local standard 
minimum of the scalar potential. 



2.3 Constraints from the dark matter relic density 

In the cMSSM, small values of mo give rise to a stable charged slepton LSP, which would 
be an unacceptable dark matter candidate. The slepton LSP problem in the cMSSM with 
small rriQ can be evaded in the cNMSSM due to the presence of the additional singlino-like 
neutralino which, in large regions of the parameter space, is the true LSP [7]. However, in 
order to be a good dark matter candidate, its relic density should comply with the WMAP 
constraint [8] 

0.094 < n^oy < 0.136 (at 2a) . (2.16) 

In the case of a singlino-like LSP, the upper bound on the relic density implies that the 
singlino-like LSP mass m^^ (= m^o) has to be close to (but somewhat below) the mass 
of the NLSP, which in the present case is always the lighter (mostly right-handed) stau 

(2-17) 

Only then can the singlino co-annihilate sufficiently rapidly with the NLSP. The con- 
dition for nearly degenerate stau and singlino masses can be obtained by replacing in the 
singlino mass squared, m^^ ~ 4fi;^s^ [1-3,9,17], the analytic approximation for the singlet 
vev as obtained from eq. 02.91) . 

1 

Ak 



-A. + y'A^ - 8m|J . (2.18) 
Noticing again that ~ Aq (< 0) and m| ~ mg, m^^ is approximately given by 

^ (^Al + \Ao\ ^Al-8ml^ - 2ml ■ (2-19) 



ml^ ~ 



An analytic approximation for the right-handed stau mass at the weak scale, obtained 
by integrating the RG equations, is given by [2] 

ml^ ml + 0.1 M^/^ . (2.20) 
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Inserting eqs. fl2.19p and (12.201) into eq. fl2.17p and using eq. fl2.1ip . one can derive the 
bound 

ml < ^Ml,^. (2.21) 

In practice, however, all approximations above (notably the neglected tr — fi mixing) tend 
to overestimate mo, and the stronger bound 

mo < (2.22) 

holds. Hence, mo must be quite small when compared to M1/2, and could well vanish. 

These analytic approximations allow to understand which "hyperplane" in the para- 
meter space Ml/2, mtQ and Aq will satisfy the WMAP constraint of eq. fl2.16p : not only 
must mo be small, but Aq is essentially determined in terms of M1/2 by eq. (I2.17p . Our 
numerical results (for mo ~ 0, see below) correspond to 

Ao ~ . (2.23) 

Finally, as in the cMSSM [22], the WMAP constraint also requires that M1/2 should not 
be too large. Since the dominant annihilation process of R-odd sparticles is now fi + ti ^ 
SM particles, the rate decreases with m^^ (which is roughly proportional to M1/2), eventually 
becoming too small for M1/2 > 2-3 TeV. 

The routine NMSSMTools [9], which includes a version of the dark matter tool 
MicrOMEGAS [10] adapted to the NMSSM, allows to scan the parameter space of the 
cNMSSM and to verify which parameters M1/2, ""^o and Aq satisfy the WMAP constraint, 
eq. (12.161) . We recall that A is another free parameter of the model, whereas tan/5 is fixed 
by the condition ms = mo at the GUT scale. In the next subsection we will discuss that 
LEP constraints impose an upper bound on A, A < 10~^, and, for illustrative purposes, 
we will fix A = 0.002 throughout the remaining part of this subsection. In any case, the 
following results are practically independent of A. 

In Fig.[T]we show the allowed range of the parameters mo and Aq as functions of Mi/2- In 
the left panel, we display the contours for the minimal and maximal values of mo. We notice 
that mo = is only compatible with the WMAP constraint for M1/2 ^ 2 TeV. Moreover, 
as Ml/2 increases, it becomes increasingly difficult to satisfy the WMAP constraint and, 
hence, the allowed range for mo decreases. For M1/2 > 3 TeV, hardly any parameter space 
survives since ^y^oh? would be too large. 

In the right panel of Fig. [H we present the values of Aq corresponding to the maximal 
and minimal values of mo for a given M1/2 (we recall that eq. fl2.16p fixes Aq in terms of 
Ml/2 and mo). Likewise, we display the values of tan/5 obtained from the requirement of 
full scalar mass unification m| ~ rn^. 

We notice that tan /? turns out to be quite large (see Ref. [23] for earlier work on the 
NMSSM at large tan /5). The origin of the large value of tan (3 can be understood as follows. 
First, an effective 5-parameter 

Eeff = Ax + KS (2.24) 

can be defined, which plays the same role as the 5-parameter of the MSSM. Then, tan/5 
is inversely proportional to ~ |-Beff|- In the regime where m| ~ ml <C ~ A^, eq. (12.180 
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Figure 1: In the left panel, maximal and minimal values of tuq allowed by the WMAP 
constraint of eq. (12.161) . as a function of M1/2 (in GeV). In the right panel, associated range 
of Aq (in GeV) and tan/3, also as a function of M1/2. In the case of tan/3 (^o), the lower 
(upper) line corresponds to mo = 0. 



gives s ~ —Ai^/2k, and thus i?cfr ~ ~ with ~ Aq and Ax determined by the RG 
equations. Finally, accidentally. Ax — \Ai^ happens to be small (much smaller than /igg) 
leading to tan/3 ^ 1. 

The lower limit on M1/2 of ~ 400 GeV follows from the lower bound on the lightest stau 
mass of ~ 100 GeV from the negative LEP searches [19]. The corresponding lower limit on 
Ml/2 derived from eq. fl2.20p seems somewhat weaker, but the stau mixing has to be taken 
into account. For M1/2 ~ 400 GeV, we observe from Fig.[T]that only mo < 20 GeV is viable. 
For larger values of M1/2, values for mo up to ~ ^Mi/2 are possible (with sizable values 
of tan/3 > 40). For values of M1/2 larger than 2 TeV, compatibility with the correct dark 
matter relic density (which requires M1/2 < 2 TeV for mo = 0) can no longer be obtained for 
mo = 0: the increasingly larger values of tan/3, as determined by the requirement of scalar 
mass unification at the GUT scale, generate a stronger mixing in the stau sector. Then, 
in order to have m^^ > m^o, a non- vanishing, albeit small, value of mo is required. For 
mo 7^ 0, compatibility with the WMAP bound allows for M1/2 up to around 3 TeV, where 
mo ~ 300 GeV and tan/3 ~ 46, with the upper bound on mo following from eq. (12.111) . 

Restricting ourselves to the phenomenologically more interesting regime M1/2 ^ 1-5 TeV, 
we present in Fig. [2] several WMAP-compatible "lines" in the [Mi/2,^0] plane for fixed 
values of mo. Actually, the allowed 2 a range for Vlh? implies that the lines displayed 
correspond to "bands" of finite (yet small) width. The parameter space lying below these 
lines is typically excluded due to the presence of a stau LSP (and to a minor extent, also a 
violation of the WMAP constraint). The upper regions delimited by each line are excluded 
due to an excessively large relic density. 
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Figure 2: Distinct allowed "lines" in the [Mi/2,^0] plane compatible with the WMAP 
constraint of eq. (12.161) . From top to bottom, the lines correspond to decreasing values of 
mo = from 20 to 140 GeV. 



2.4 Constraints on A from LEP and dark matter 

Remarkably, LEP constraints on the SM-like Higgs boson mass turn out to be satisfied due 
to the relatively large stop masses and trilinear coupling At for all parameter ranges shown 
in Figs. [1] and O However, they lead to upper bounds on the NMSSM specific parameter A. 
For the large values of tan /3 obtained in this scenario, a large value of A does not lead to 
an increase of the SM-like Higgs mass. On the contrary, increasing A simply increases the 
mixing of the singlet-like CP-even scalar with doublet-like CP-even scalars. If the singlet- 
like CP-even scalar mass is larger, the SM-like Higgs mass decreases with increasing A and 
can fall below the LEP bound. If the singlet-like CP-even scalar mass is below the SM-like 
Higgs mass limit, i.e. < 114 GeV, its coupling to the Z-boson, which is proportional to A, 
must be sufficiently small, equally implying an upper bound on A. 

In Fig. [3], we show the corresponding upper limits on A for the case mo = 0. The 
constraint is particularly strong in the "cross-over" region near M1/2 ~ 660 GeV (see Fig. [5] 
below), where relatively small values of A can generate a large mixing angle. We see that 
A < 10~^ is required for all values of M1/2; hence, the parameter A will have practically no 
effect on the remaining spectrum except for the singlet-like Higgs masses. 

On the other hand, a rough lower bound on A can be derived from the efficiency of 
singlino-like LSP annihilation in the early universe. For the dilution of the LSP, two pro- 
cesses are relevant in the limit of small A. The dominant annihilation process of R-odd 
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Figure 3: The upper bound on the parameter A as a function of M1/2 (in GeV), as obtained 
from LEP constraints on the NMSSM Higgs sector; for simphcity, we set mo = 0. 



particles is NLSP + NLSP X, the NLSP being the hghtest stau. The reaction rate for 
this process is given by (riNLSp)^cr, where tinlsp is the corresponding time and temperature 
dependent abundance, and a is the thermally averaged cross section. 

The other relevant process leading to the observed relic density is LSP+X NLSP+X' 
(and its inverse), where X and X' are practically massless quarks and leptons - this process 
helps to maintain the LSP and the NLSP in thermal equilibrium. Its reaction rate can be 
written as nLSP nx (j', and depends on the small, but non- vanishing, non-singlet component 
of the singlino-like LSP of 0{X/g2), where (72 is the SU(2) coupling constant. Consequently, 
the cross section a' is ~ A^/f^l x a, and is correspondingly suppressed for small A. 

Nevertheless, the process LSP + X —>■ NLSP + X' is typically faster than the annihilation 
process NLSP + NLSP —>■ X, since near the freeze-out temperature, the abundance nx 
of quarks and leptons is ~ 10^ larger than the abundances of the LSP and NLSP (for 
mLSP ~ ■'^NLSp) [24]. This allows to dilute the LSP density as fast as the NLSP density. 

Only for very small A, the reaction rate of the process LSP + X — >■ NLSP -|- X' can 
become smaller than the one for NLSP + NLSP —* X; then the LSP will no longer be in 
thermal equilibrium with the NLSP near the freeze-out temperature, but can be considered 
as decoupled, implying an excessively large relic density. According to the discussion above, 
this would happen for \^ / ^ 10~^ or A < 10~^. Clearly, these are rough estimates, which 
would merit more detailed investigations. In the following, we employ values of A sufficiently 
above 10~^ such that the hypothesis of thermal equilibrium between the LSP and the NLSP 
near the relevant temperature can be considered as satisfied. 

We remark that, even if we hypothetically allowed for other contributions to the dark 
matter relic density, this would not affect the lower bound for A estimated above, nor the 
previous discussion on the allowed parameter space (i.e., the derived bounds on M1/2, 
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etc.), since all bounds originate from the upper WMAP limit on the relic density. 



2.5 Constraints from flavor physics 

As previously mentioned, we have checked that constraints from 5-meson physics [25] are 
satisfied. More precisely, agreement within 2 a between the following observables and their 
theoretical values is verified by the NMSSMTools routine [9]: the decay branching ratios 
BR{B Xs-f), BR(i?+ T^i^r), BR{Bs Ai+/i") and the mass differences AMq, q = d,s. 
It turns out that all regions in parameter space investigated before (i.e. consistent with 
WMAP and collider constraints) are also allowed by constraints from i?-physics. 
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Figure 4: 5aSUSY 

function of Mxji (in GeV) for mo = 0. The vertical bars denote the 
theoretical error, while darker (lighter) regions correspond to a 1 cr (2 cr) deviation from the 



central value of a^f^ 



,SM. 



data taken from Ref. [27]. 



More important constraints arise from requiring that SUSY accounts for the ~ 3 a 
deviation of the observed anomalous magnetic moment of the muon with respect to the SM 
expectatioifl, 5aJ"^^ = (27.7 ±9.3) x lO^^o [15]. The cNMSSM has been analysed in this 
respect in Ref. [27], the analysis being conducted for simplicity for mo = 0, but this result 



is practically independent of mo. In Fig. HJ we show ba 



SUSY 



[g — 2)^^^^ as a function of 



Ml/2, depicting the 1 a and 2 a bands as well. 

We see that the constraint from would confine the allowed range of M1/2 to M1/2 ^ 
1 TeV at the 2 a level, and to 400 GeV < M1/2 < 600 - 700 GeV (where the sparticle 
spectrum is not too heavy) at the 1 a level. In fact, the present experimental value could 
be matched to arbitrarily high precision, and a more precise measurement of would 
eventually lead to a prediction of M1/2 in the cNMSSM. In any case, in view of the desired 
value for (5a^^^^, the region M1/2 ^ 1 TeV is preferred. 



recent measurement of the hadronic cross section e^e — > hadrons using radiative return by the 
BaBar collaboration indicates that this discrepancy might be smaller than presently thought [26]. 



11 



3 Higgs and sparticle spectra 



We now proceed to analyse in some detail the Higgs and sparticle spectra obtained in the 
allowed regions of the cNMSSM parameter space, after imposing the bounds from Higgs 
and sparticle searches at LEP, the 5-meson constraints and the requirement of a correct 
cosmological relic density as measured by WMAP for the neutralino LSP. 

3.1 The Higgs spectrum 

The Higgs sector of the (c)NMSSM contains three neutral CP-even Higgs states /i^ {i = 
1, 2, 3), two neutral CP-odd Higgs states a° (i = 1, 2), and the charged Higgs states /i^. 
Their masses depend essentially on the gaugino mass parameter Mi/2- 

As in the MSSM in the decoupling regime (see Ref. [28] for a review), the heaviest CP- 
even, CP-odd and charged Higgs states form a practically degenerate SU(2) multiplet with 
a common mass above 500 GeV. The mostly SM-like CP-even state has a mass increasing 
slightly with M1/2 from 115 GeV up to ~ 120 GeV. This mass range is only slightly above 
the lower limit of 114.4 GeV on the SM Higgs boson mass, and is compatible with the Higgs 
mass range favored by electroweak precision data as recently obtained from a global fit (in 
which the central mass value is 116 GeV) [29]. 

The third CP-even state has a dominant singlet component; it is the only Higgs state 
whose mass depends - apart from M1/2 - on mo and, to some extent, on A: for small M1/2 
it is lighter than the SM-like Higgs boson, escaping LEP constraints due to the very small 
coupling to the Z boson. For increasing values of M1/2, its mass increases until it becomes 
nearly degenerate with the SM-like CP-even Higgs state: in this region of parameter space, 
which will be subsequently denoted as the "cross-over" region, the singlet-like and SM-like 
Higgs states strongly mix. For still larger values of M1/2, the mass of the singlet-like state 
exceeds the one of the SM-like state. 
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Figure 5: The Higgs masses as a function of M1/2 (in GeV). In the left panel we set mo = 0, 
while in the right panel mo is given by its maximal value, m™'^^(Mi/2). From below, the 
displayed lines correspond to the states (blue/dotted), /i^ (pink/dashed), a° (full/black) 
and (full/red) which is degenerate with the and states. 
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This "cross-over" phenomenon is visible in Fig. O where we display the masses of the 
neutral CP-even, CP-odd, and charged Higgs bosons as a function of the parameter Mi/2- 
On the left-hand side, we take tuq = 0, while on the right-hand side we assume mo = 
m™'^^(Mi/2), as given by the upper line in the left panel of Fig. [H (With the exception of 
the CP-odd state, the Higgs mass spectrum is somewhat heavier in the case mo = 0.) 

In Fig. [5] we have set A = 0.002; then, for mo = corresponding to the left panel, 
the "cross-over" phenomenon occurs at M1/2 ~ 660 GeV: for M1/2 ^ 660 GeV, the light- 
est CP-even state is singlet-like, whereas the lightest CP-even state is SM-like for 
Ml/2 ^ 660 GeV. In the case mo = m^'^^{Mi/2) displayed in the right panel of Fig. [5l the 
corresponding "cross-over" occurs at M1/2 ~ 1100 GeV. Intermediate regimes for mo (with 
mo 7^ 0) will imply a "cross-over" that takes place for 660 < M1/2 < 1100 GeV. 

For the present value of A = 0.002, the CP-even singlet mass can be as small as 60 GeV 
for mo = m™'^^(Mi/2) (in the right panel). For larger values of A, close to the upper bound 
shown in Fig. [3l even lower m/^o values can be obtained, as will be discussed later. 

The lighter singlet-like CP-odd scalar has a mass above ~ 120 GeV (increasing with 
M1/2); the heaviest CP-even and CP-odd scalars /ig and are practically degenerate in 
mass with the charged Higgs boson h^, with masses above ~ 520 GeV. 

3.2 A possible explanation of the excess in Higgs searches at LEP 

In addition to the 1.7 a signal for a SM-like CP-even Higgs particle with a mass close to 
115 GeV, the combined results on Higgs searches of the four LEP experiments via the 
Higgs-strahlung process e~^e~ —>■ hZ, followed by the Higgs decay h bb, show a 2.3 a 
excess of events corresponding to a Higgs mass around 98 GeV [12,16]. The number of 
excess events amounts to about 10% of those expected for a SM Higgs boson with the 
same mass. It can be explained either 

i) by a reduced coupling of a candidate Higgs boson to the SM gauge bosons, = 

9hzz/ghSMzz ~ O{V0A), 
a) by a reduced branching ratio of a candidate Higgs boson into bb final states. 

In the unconstrained MSSM, the two excesses can be explained [30] by the presence of 
a SM-like Higgs boson with a mass of ~ 115 GeV (the heavier CP-even H state), while 
the lighter CP-even state h has a mass close to ~ 98 GeV and reduced couplings to the Z 
boson. 

Within the context of the unconstrained NMSSM, the explanation ii) of the excess for 
Higgs masses around 98 GeV has been proposed [31]. In this case, one can have a CP-even 
Higgs boson with a corresponding mass and SM-like ZZh couplings, i.e. ~ C'(l), but 
a reduced branching ratio into bb final states due to a dominant decay into pairs of very 
light CP-odd bosons h\ — > aja^ with m^o < 2mh. Then, the light pseudoscalar a\ can only 
decay into r'^r" and eventually light quark and gluon pairs, rendering it compatible with 
corresponding searches at LEP [12]. 

In the cNMSSM, the parameter space somewhat below the cross-over regions contains 
neutral Higgs scalars with masses ~ 100 GeV and couplings ~ (9(-\/0.1). As an example, 
for 560 GeV < M1/2 < 575 GeV (with mo = and A = 0.005), one finds for the masses of 
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the two lighter CP-even Higgs states 97 GeV < rrif^o < 101 GeV and rrif^o ^117 GeV. In 
this case, the reduced couphng of the hghtest scalar Higgs boson to SM gauge bosons would 
lie in the range 0.28 < |C^o| < 0.33 for so that the cNMSSM could indeed account for 
the observed 2.3 a excess around rrih ~ 98 GeV at LEP. In addition, since the mass of the 
nearly SM-like /ig state is m^o ^ 117 GeV with ~ 0.9|, and in view of the error of a 
few GeV in the determination of the radiative corrections to the Higgs masses (expected to 
be, as in the MSSM, of the order of 3 GeV, see e.g. Ref. [32]), the 1.7 a excess at a Higgs 
mass ^115 GeV could be explained as well. 

Thus, in the cNMSSM, the region in parameter space corresponding to small M1/2 can 
describe not only the deviation of the {g — 2)^ from the SM expectation, but both excesses 
of Higgs-like events at LEP as well. 

3.3 The sparticle spectrum 

Let us now turn to the sparticle spectrum, starting with the neutralino and slepton mass 
spectra shown in Fig. [61 As for the Higgs bosons, we display the case mo = in the 
left-hand panel, and rriQ = m™'^^(Mi/2) in the right-hand panel. 

The two nearly degenerate sets of lower lines in both panels correspond to the masses of 
the Xi singlino-like LSP (blue/dotted) and the lighter stau Ti NLSP (red/full). The mass 
difference between these two states is smaller than ~ 8 GeV, as required in order to obtain 
a cosmological relic density for the singlino Xi compatible with WMAP. 

The pattern for the masses of the charginos and the heavier neutralinos (blue/dotted 
lines) follows the one of the MSSM, once the proper relabeling of the states is made. 
Since the low-energy value of the higgsino mass parameter fi^fi is generally quite large, 
A^eff ^ M2, the heavier neutralino states xt ^-nd xt ^-^e higgsino-like with masses ~ fies- 
The states X2 X3 ^-^e, respectively, bino and wino-like with masses m^o ^ 2m^o ^ M2 




Figure 6: Neutralino (blue/dotted lines), selectron (black/full lines) and stau (red/full lines) 
masses as a function of M1/2 (in GeV); on the left-hand side we set mo = 0, while on the 
right-hand side we set mo = m™'^^(Mi/2)- In both panels the states are ordered in mass 
as m^o < mf J < mg^ < m^o < mf^ < m^^ < m^o < m^o ^ . The charginos xf and x^ are 
degenerate in mass with, respectively, Xs and Xi5- 
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Figure 7: Squark and gluino masses as a function of M1/2 within the range obtained by 
taking mo = or mo = m™'^^(Mi/2)- The states are mass ordered as m^^ < m^^ < ''^wjj < 

(with M2 ~ 0.75Mi/2)- The charginos xf ^ind X2' nearly degenerate in mass with, 
respectively, the wino-like X3 ^^id the higgsino-like Xa 5 states. 

For completeness, we have also depicted the masses of the left and right handed selec- 
trons, ei and (black/full lines), which are degenerate in mass with, respectively, the 
smuons fli and fiR. The right-handed states are always lighter than the left-handed ones 
and the mass pattern is such that mg^ < m^o < mg^. Of course, all these states are heavier 
than the NLSP fi (in fact, is almost degenerate with 7^2 ). Due to SU(2) symmetry, the 
sneutrinos have approximately the same masses as the left-handed sleptons. 

In Fig. [TJ we display the masses of the gluino and of up-type squarks as a function of 
Ml/2. Here, the variation of the parameter mo in the range < mo ^ m™'^^(Mi/2) results 
only in a hardly visible "width" of the lines. The masses of the c-quarks are degenerate 
with those of the u-squarks and have not been displayed in the figure. 

We first note that the gluino is always heavier than the squarks, a consequence of the 
small ratio mo/Mi/2; this feature has important consequences as will be discussed later. The 
left- and right-handed up-type squarks are almost degenerate (the mass difference being 
less than 5%), and also nearly degenerate with the first and second generation down- type 
squarks, which are not shown in Fig. [71 The top squarks are lighter than the u-squarks, the 
mass of the lighter stop ti being ^ 20% smaller. Thus, ti is the lightest strongly interacting 
particle. The bottom squarks have masses somewhat below the heavier stop ^2- 

Actually, since the singlet sector practically decouples, most of the spectrum could also 
be obtained in the cMSSM, provided the present values of M1/2, mo and Aq are used. We 
have explicitly verified that the obtained spectra of the non-singlet states for the two points 
of the cNMSSM parameters discussed in the next subsection agree with those obtained in 
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the MSSM with the program Suspect [33] (up to small differences due to the treatment of 
the higher order radiative corrections). 

3.4 Examples of spectra 

The most relevant features of the sparticle and Higgs spectrum can be represented by 
two points PI and P2, for which the distinct cNMSSM inputs, as well as the spectra, are 
summarized in the left part of Table [H These two points illustrate the low and intermediate 
Ml/2 regime in the cNMSSM. 

In the right part of Table [H we display the resulting spectrum for the Higgs, LSP and 
NLSP states for rriQ ^ 0. While the input M1/2 is the same, the input values of Aq and 
tan P have been slightly adjusted in order to obtain acceptable values for the dark matter 
density and the unification of the singlet mass. In the points PI and PI', the CP-even 
state is singlet-like, with a somewhat lower value of mj^o for PI'. The state is SM-like 
for P2, while hi and h2 have a very similar non-singlet components in P2' implying similar 
couplings to gauge bosons and quarks. 

4 Prospects for collider searches 
4.1 Sparticle and Higgs decays 

The most interesting aspects of the spectrum of the cNMSSM - as compared to the cMSSM 
- are the presence of a singlino-like LSP with a mass just below the one of the stau NLSP, 
and the fact that all squarks are lighter than the gluino. The singlino-like LSP with its 
small coupling to all other sparticles will strongly modify the sparticle decay chains, since 
now all sparticles will decay via the stau NLSP. Also, squark and gluino decay chains will 
differ from most MSSM-like scenarios. These two features will have important consequences 
for sparticle searches at the LHC 

4.1.1 Sparticle decay branching ratios 

The branching ratios for the gluinos, squarks, sleptons and charginos/neutralinos are shown 
in Table [2] for the points PI and P2 with mo = (the situation being similar in the primed 
points with mo 7^ 0). They have been obtained by applying the program SUSYhit [34], 
which calculates the decay widths and the branching ratios of the Higgs and SUSY particles 
of the MSSM, to the NMSSM with a practically decoupled singlet sector. A few comments 
are in order. 

- As the gluino g is heavier than all squarks, it can decay via two-body decays into all 
quark-squark pairs. The branching ratio into tti final states is somewhat larger (~ 20%) 
as a consequence of the larger phase space due to the lighter ti states. 

- All squarks (including the stops) decay into neutralinos or charginos plus the corre- 
sponding quark. For right-handed squarks qn of the first two generations, the branching 
ratio of the decay qn ^ X2Q nearly 100%, the X2 state being dominantly bino-like. In 
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Table 1: Input parameters and low-energy spectra for four points of the cNMSSM with two 
distinct M1/2 regimes. On the left are Points PI and P2 with mo = 0, M1/2 = 500 GeV 
and 1 TeV. On the right PI' and P2' with mo 7^ 0, for which we only display those masses 
whose values vary with niQ. In all cases, we have set A = 0.002. 

the case of the left-handed states, the branching ratios are ~ | and ~ | for the decays 
into the neutral qi — > qxt charged qi — > q'Xi wino states, respectively. 

- Regarding the decays of the electroweak gauginos, the preferred decay channel of the 
state X2 is X2 ~^ ^1'^ which has a branching ratio of ~ 90% (for PI) as a result of the 
more favorable phase space; the remaining ~ 10% are the decays xt ~^ ^1 where / = 
or /i^. The wino-like Xg and xt decay ~ 50% into first/second generation slepton-|-lepton 
states and ~ 50% into third generation fr, VrVr states; the reason for the breaking of lepton 
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Table 2: Dominant decay modes of the squark, slepton, gluino, neutralino and chargino 
states for the two points PI and P2, for which the spectrum is given in Table [1] (and where 
A = 0.002). They have been obtained using SUSYhit [34]; q and / denote first and second 
generation quarks and leptons, respectively. 

universality is again the more favorable phase space. 

- Finally, while the left-handed first /second generation sleptons II decay to 100% into 
1X2 final states, the right-handed sleptons Ir essentially decay via the three-body channel 
k ~^ l^iT'i this decay mode has also been discussed in Ref. [35], albeit in a different 
context. The branching ratio for the two-body decay mode — ^ Xi^ is well below the 
percent level; the decay into the bino X2 ^"^^ ^ lepton is forbidden by phase space. The 
reason for the dominance of the three-body decay is that the two-body decay can occur 
only via the bino-component of 0{X) of the mostly singlino-like Xi? and hence is extremely 
small, even for the maximally possible values of A ~ 0.01 shown in Fig. [31 On the other 
hand, the three-body decay occurs through the virtual exchange of the bino-like whose 
virtuality is not very large as the Ir and X2 masses are comparable. 

Hence, practically all sparticle decay chains contain the fi NLSP. The fi life time can be 
very large, but it will finally decay into the singlino-like LSP and a tau lepton, fi — > Xi 

4.1.2 Displaced vertices 

For very small A, the couplings between the fi, x\ ^"^^ states might be sufficiently small, 
resulting into a stau track of O(fewmm) that might be visible [36,37]. Hence, displaced 
vertices at high-energy colliders such as the LHC from long-lived staus could be a "smoking 
gun" signature of the cNMSSM. Here we present some details of the computation and the 
resulting possible track lengths. 
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Following Ref. [38], the partial width of the stau decay into the lightest neutralino and 
tau lepton can be written as 



Tih^XiT) = ^[{al,+bl,){ml-ml-mlo)-ianbumrm^o], (4.25) 



Tl 



where p is the phase-space function, y, z) = + y"^ + — 2xy — 2xz — 2yz, and an, 
bii are the Xi tr,l couplings, defined as 

O'li = " g2V2 tan 6^Jv Nu sin Or — h^- N13 cos Or , 

611 =^(iVi2 + tan^wA^ii) cos^^ - K N13 sinOr- (4.26) 
v2 

hr and Or denote the tau Yukawa coupling and stau mixing angle, (?2 the SU(2) coupling 
constant and A'^ij denote the composition of the lightest neutralino, Xi = ^iiB + N12W + 
NisH^ + Ni4^H^ + N15S. In the regime of small A, the bino (A^n), wino (A''i2) and higgsino 
(iVi3 and A'^14) components of the mostly singlino-like neutralino Xi ^-^^^ all proportional 
to A [36]. Hence it is convenient to introduce nearly A-independent coefficients a and (3 
defined as af-^ + bf-^ = \^ a and 2aii6ii = A^/5. 

In the relevant limits Am = m^^ — m^o <^ m^^ ~ m^o and rrir <^ ~ ""^x?' 
expression eq. fl4.25p for the stau decay width can be simplified as [36] 



A^Am^ 



m. 



2 



m - X?r) ^ A^ ^ ^ («Am - /5m.) , (4.27) 

which summarizes the essential dependence of the stau decay width on A and Am. 

The coefficients a and /? still depend somewhat on M1/2 and mo, and their numerical 
values smoothly decrease with M1/2; for 400 GeV < M1/2 ^ 1500 GeV, one has 0.01 
> a ~ /5 > 0.0001. The mass splitting Am does not only depend on M1/2 and mo: the 
allowed 2 a range for the relic density in eq. f l2.16p allows Am to vary within a small window 
at fixed Ml/2 and mo, but slightly fiuctuating Aq. In the left panel of Fig. [HI we show the 
corresponding ranges for Am. Note that for Am < m. the decay of the fi has to proceed 
via a virtual r, implying a tiny partial width. 

On the right panel of Fig. [Hlwe plot the reduced stau length of fiight, ll'f' = hc/V{ji — > 
Xit), as a function of M1/2 for A = 10~^. For other values of A, the reduced length of fiight 
can be obtained by rescaling Vf^ by a factor (10~'^/A)^. For M1/2 > 1200 GeV and mo = 0, 
the lifetime of the stau can be extremely large as a consequence of the NLSP-LSP mass 
difference approaching the rrtr threshold, which corresponds to the vertical dotted line that 
extrapolates the upper-most curve in the right panel of Fig. [H 

The stau length of flight in the laboratory frame is given by 



W = /^f V^/5|,/(1-/5|J, (4.28) 



where (3f^ = Vf^/c is the ri velocity. A realistic estimate of If^ requires the knowledge of (3f^ 
and hence of the production processes of the lightest stau. 
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Figure 8: On the left panel, the range of neutralino-stau mass differences Am which are 
allowed by the 2 a error bars of the WMAP constraint, as a function of Mi/2- The lower 
lines correspond to rriQ = 0, the upper ones to its maximal possible value. In each pair, full 
(red) lines denote the maximal Am - associated ^^oh'^lma.x, while dotted (blue) correspond 
to f2^o/i^|inin- The horizontal line denotes m^. On the right panel, the maximal and minimal 
"reduced" stau lengths of flight, If'^'^ (in mm), as a function of M1/2, for A = 10~^, for the 
neutralino-stau mass differences given in the left panel. 



A hint on realistic values for can be obtained from the gauge-mediated SUSY break- 
ing ATLAS benchmark point GMSB5, where it is advocated that over 99% of the staus 
(which will decay into a gravitino LSP) have ^0.7 [39]. 

In Fig. [9l we display the decay length If^ as a function of M1/2 for the choice mo = 0, 
taking an intermediate value of Am within the range presented on the left panel of Fig. [8l 
We consider two values of and two distinct regimes for A. As mentioned above, even 
for A close to its upper limit (left panel of Fig. [9]), visible lengths of flight If^, of (9(mm), 
are possible for the lightest cNMSSM stau for large M1/2 ~ 1-4 TeV. 

On the right panel of Fig. [9] we consider a smaller value A ~ 10~^, and focus on the 
possibilities within the range of M1/2 < 1 TeV favored by the data on [g — 2)^, where one 
can still obtain a stau length of flight as large as a few centimeters. For larger values of 
Ml/2 and A < 10"^, If^ can become as large as C(lOcm) and even 0(lm). 

Long-lived stau NLSPs are also possible in MSSM-like models, if the gravitino is the 
true LSP. However, at least if m3/2 is not much smaller than MgusY, the stau length of flight 
will always exceed the size of the detector (see [40] for recent analyses of collider signatures 
of such scenarios). For gauge-mediated SUSY breaking models, where m3/2 ^ Mgusy, 
implications of long-lived stau NLSPs for SUSY searches at high energy colliders such as 
the LHC have been discussed for example Ref. [39]. For the NMSSM, no such studies exist 
at present. 

In any case, although visibly displaced vertices from long-lived stau decays are not an 
unavoidable prediction within the cNMSSM, they are a very interesting possibility of the 
parameter space associated to small A and/or Am. 
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Figure 9: The stau length of flight /^^ (in mm), as a function of M1/2 (in GeV), for niQ = 
and A = 10"'^ (left panel) and A = 10^"^ (right panel). In both panels, the upper line 
corresponds to /^^^ = 0.95, while for the lower one we have chosen = 0.7. Note the 
different scales of the Mi/2-axis. 



4.1.3 Higgs decays 

In the cNMSSM, the decay pattern of the non singlet Higgs particles will follow that of the 
MSSM in the decoupling regime [28]. 

The lightest non-singlet like CP-even Higgs boson is SM-like with a mass between 115 
and 120 GeV, implying that it will dominantly decay into bb pairs with a branching ratio 
larger than ^ 70%, followed by decays into cc, t~^t~ and gg pairs, with branching ratios 
of the order of 5%. The branching ratio for the decay into a pair of W bosons (one of 
them being virtual) is less than 10%. The branching ratio for the interesting decay into 
two photons will be at the level of BR(/i^^ ^7 7) ~ 2 x 10~^. 

The heavier non-singlet Higgs particles /13, 03 (and h^) will mainly decay into bb and 
T~^T~ {tb and tu) pairs. The branching ratios will be 90% and 10% for the hadronic and 
leptonic decay modes, respectively. Due to the strong enhancement of the b and r Higgs 
couplings as a result of the large value of tan /3 > 30, the Higgs decays into SUSY particles, 
as phase space allowed decays into pairs of the lighter neutralinos and sleptons, will be 
strongly suppressed. 

The decays of the singlet-like neutral Higgs bosons (which are unlikely to be produced as 
will be discussed later) are induced through their (typically small) non-singlet components, 
implying branching fractions corresponding to a SM-like Higgs. 

Higgs-to-Higgs decays, which are a very peculiar possibility within the unconstrained 
NMSSM and which have been discussed in great detail in Ref. [5], are generally absent in 
most of the parameter space of the cNMSSM; only for the largest possible values of A, the 
lightest CP-even singlet-like Higgs mass can be sufficiently small and its coupling to still 
sufficiently large as to give rise to a non- negligible BR(/i2 — > hi hi). An illustrative example 
is given by the following choice of input parameters 

Ml/2 = 500 GeV, mo = 46.5 GeV, Aq = -142.65 GeV, tan/? = 31.72, A = 0.0097, 

(4.29) 
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where one obtains m/jO ^ 50 GeV and m^jO ~ 118 GeV. The dominant /ij branching fractions 
are then 



BR(/i° ^bb) = 64.1% , BR(/i° ^ /i° hi) = 12.1% . (4.30) 



The state will mostly decay into bb (90%) and t^t (10%) final states. 
4.2 Sparticle and Higgs production at the LHC 

Regarding sparticle production at the LHC, we summarize in Table [3] the cross sections for 
the most relevant production channels for the points PI and P2 defined in Tabled] These 
have been obtained by using the Fortran code Prospino [41], which calculates the cross 
sections for sparticle production at hadron colliders to next-to-leading order in perturbation 
theory, for the non-singlet sector of the cNMSSM. 

The largest cross sections are expected for the strongly interacting production processes 
pp — > qg, qq and qq*, where q denotes all squarks but the stop. One obtains agg ~ 0.67 pb, 
agg ~ 0.44 pb and agg* ^ 0.22 pb for the points PI (and PI') with M1/2 = 500 GeV, 
while for the point P2 (likewise for P2') with M1/2 = 1000 GeV, only agg ^ 0.01 pb, 
agg = 0.004 pb and agg* = 0.002 pb are expected. Assuming the LHC high luminosity 
option of £ = 100 fb~^, one obtains 10^ to 10^ events for scenario PI, but only 10^ to 
10^ events for scenario P2 (before efficiency cuts are applied). The yield for gluino pair 
production, qq/ gg —>■ gg, is much smaller as a result of the larger gluino mass. 

As discussed in Section 4.1, all gluinos will decay into squarks. Squarks decay into 
quarks and the wino-like neutralinos or charginos, which then cascade mostly into the 
lighter f NLSP. Hence all sparticle decay chains contain the fi NLSP, which will finally 
decay into the singlino-like LSP, fi — > Xi '^^ possibly leading to a displaced vertex. 

The simplest possible squark decay cascades at the LHC - originating mostly from 
right-handed squarks of the first generation q^ - are thus of the form 



where, according to Table [2], the first case occurs ~ 88% for point PI and ~ 74% for point 
P2. The second case occurs only ~ 12% and ~ 25% for PI and P2, respectively. These 
cascades typically lead to events with 3 jets per q^ (one hard quark and two tau jets, 
also potentially hard, depending on the momentum of the decaying X2 neutralino), but the 
leptonic decays of the r can possibly also be used. 

As can be seen from eq. (14.311) . right-handed squark cascade decays containing a lepton in 
the final state are clearly sub-dominant. In the case of left-handed squarks, the decays would 
lead to more e,/i lepton final states, since qi mainly decays into qxl and qx^- The wino- 
like neutralino X3 and chargino Xi" decay either directly to fi and z>r (the latter decaying 
predominantly into fiW) or to Z^, z/;. These then decay into 1x2 1^X2^ respectively. The 
bino-like X2 would dominantly lead to fir and, to a minor extent, to IrI Utit final 
states. Therefore, one has a non- negligible probability for one, two, three and even four e, /i 
leptons in the final states. For instance, the branching ratios into the four lepton topology 
would be ~ 1% and ~ 2% for, respectively, points PI and P2. 




\ qlRl qlfiT 



qrrxl (PI : 88%; P2 : 74%) 

qrrxl + i (PI : 12%; P2 : 25%) 



(4.31) 
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a (pb) 


PI 


P2 


99 


9.5 X 10-2 


2.14 X 10-^ 


gq 


0.668 


4.28 X 10-3 


QQ 


0.436 


9.21 X 10-3 


qq* 


0.221 


1.64 X 10-3 


il il 


3.69 X 10-2 


2.63 X 10-^ 


IlII 


3.4 X 10"3 


1.62 X 10-3 


Ir l*R 


1.17 X 10-2 


8.87 X 10-^ 


~ ~ * 


3.58 X 10-3 


1.53 X 10-^ 


riff 


4.8 X 10-2 


3.46 X 10-3 




1.1 X 10-3 


6.22 X 10-5 


xlxl 


1.73 X 10-^ 


8.67 X 10-6 


xlxt 


5.37 X 10-^ 


6.53 X 10-5 


xUl 


1.79 X 10-3 


5.74 X 10-5 


xlxt 


6.51 X 10-2 


7.49 X 10-3 


XiXi 


3.53 X 10-2 


1.17 X 10-3 



Table 3: Production cross sections (in pb) for strongly (upper part) and weakly (lower part) 
interacting sparticles at the LHC in the points PI and P2 defined in Table [H as obtained 
with the program Prospino [41]. Here q denotes all squarks but the stop, and Z = // or e. 

For chargino and neutralino or mixed pair production, the cross sections are also shown 
in the lower part of Table (3] they are sizable enough only for point PI, where the phase 
space is not too penalizing. In the case of XiXi ^"^^ XzXii they lead to a few thousand 
events for an integrated luminosity of £ = 100 fb-^, before cuts are applied. However, since 
charginos and neutralinos also dominantly cascade into f lepton final states, only a very 
small fraction leads to the nice signature of multi e, lepton events. 

The cross sections for the Drell-Yan production of slepton pairs is of the same order 
as the one of charginos and neutralinos, the highest one being pp fif^ pair production 
(with 0" ~ 50 fb for PI, where m^^ 130 GeV) as a result of the more favorable phase 
spac^. 

Since cascade decays leading to e, /i lepton final states will be generally rare in the 
cNMSSM, the measurements of the sparticle masses from kinematical endpoints and, thus, 
the determination of some of the soft-SUSY breaking parameters, cannot be performed 

^In fact, at the Tevatron, the only process which might have a significant production cross section is 
the pair production of the hghter stau: one would have a{pp fiff ) « 15 fb for rrif^ « 100 GeV. The 
final state would then essentially consist of missing energy and two tau final states, the latter potentially 
decaying into muons with large impact parameters, since the fi lifetime can be e.g. of the order of 20 ps 
for nif^ w 100 GeV and A w 10-^. Such events would share some (but definitely not all) of the peculiarities 
of the displaced multi-muon events recently reported by (part of) the CDF collaboration [42]. However, in 
our case only a handful of events would have been produced and the r multiplicity would be far too small. 
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with the same level of accuracy as it would be the case for the MSSM. Note also that, in 
most cases, one of the leptons originates from the three-body Ir /fir decay. 

Hence, a SUSY signal can certainly be observed at the LHC by looking, for instance, 
for final states with hard jets plus a large amount of missing energy. However, precision 
measurements through the endpoints of decay spectra will be more complicated to perform 
in the case of r final states which are overwhelming in the cNMSSM. 

Finally, let us make a few comments on Higgs detection in this scenario. As previously 
discussed, in most cases the cNMSSM Higgs sector reduces to the one of the MSSM in the 
decoupling regime. The SM-like CP-even Higgs boson with its mass in the 115-120 GeV 
range can be discovered first at the Tevatron (if enough integrated luminosity is collected) 
in the Higgs-strahlung process pp VT+Higgs leading to lubb final states. At the LHC, 
the most relevant production and decay processes will be the gluon-gluon and vector boson 
fusion, gg Higgs and qq gg+Higgs, with the Higgs boson decaying into two photons 
(and possibly t~^t~ final states in the vector boson fusion process). 

We recall the phenomenon of the "cross-over" for certain values of M1/2 and mo, which 
corresponds to two nearly degenerate CP-even Higgs bosons sharing their couplings to the 
SM gauge bosons, quarks and leptons. This is the case illustrated by the last point P2' 
of Table [T], where the lightest Higgs bosons and /ig have very similar singlet / doublet 
components ~ 0.7, i.e. similar production cross sections and branching fractions, and a 
mass splitting of about 10 GeV. In principle, both states could be observecif]. For smaller 
values of A, the mass splitting in the cross-over region can be considerably smaller, less 
than 1 GeV for A < 10"^. Then, the sum of both Higgs bosons gives rise to a single SM-like 
Higgs boson peak, which can be quite difficult to resolve at the LHC. 

The heavier non-singlet Higgs particles have significant couplings to SM down-type 
fermions in view of the large values of tan/3 > 30. They can be observed in associated 
production with bb pairs and decays into r+r^ pairs for the neutral particles, and associated 
production with tb pairs with decays into ru in the case of charged Higgs bosons. However, 
this is only possible for low values of M1/2 < 700 GeV, which lead to light enough Higgs 
particles with sufficient production cross sections. In any case, a very high luminosity is 
required to observe these particles. 

Apart from the "cross-over" region, the singlet-like Higgs states are generally inaccessible 
as they couple too weakly to SM fermions and gauge bosons. Only if the mass of a singlet- 
like CP-even Higgs boson is sufficiently small and A sufficiently large, as is the case for 
the point given in eq. (14.291) . the branching ratio /ig — >■ h^h^ of the then SM-like /ig can be 
sizable (leading to a slightly reduced /ij —>■ 77 branching ratio). Whether this production 
of hi leading to difficult 46, 2r2b and 4r final state topologies could be detected [5], needs 
also a dedicated detailed investigation. 

■^The situation here is similar to the MSSM intense couphng regime discussed in Ref. [43] in which all 
neutral MSSM bosons are close in mass; however, in the cNMSSM, the couplings of the lighter CP-even 
states to down-type fermions are not particularly enhanced. 
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4.3 Searches at the ILC 



Due to the relatively heavy Higgs and sparticle spectrum, a multi-TeV e+e~ collider would 
be required to produce all the states of the constrained NMSSM. At a 500 GeV ILC, 
only the two lighter neutralinos and the right-handed sleptons can be produced even for 
^1/2 ^ 500 GeV. The production of wino-like charginos and neutralinos as well as the 
left-handed sleptons would need a larger center of mass (cm.) energy of ~ 1 TeV. The 
detection of the various SUSY states would be straightforward in the clean environment 
of the ILC. In particular, a fi nearly degenerate with the LSP can be detected, as shown 
in detailed simulations for somewhat similar MSSM scenarios [44,45]. The masses of the 
SUSY particles could be accurately determined, at least through threshold scans, as has 
been discussed in detail in Ref. [45], and a clear distinction between the cNMSSM and other 
scenarios could be made. 

As far as the Higgs sector is concerned, only the lighter CP-even and CP-odd Higgs 
states would be kinematically accessible at an ILC with a cm. energy less than 1 TeV. 
The SM-like Higgs particle can be easily detected and its properties probed in detail, the 
Higgs mass range around ~ 120 GeV being the ideal one for a 500 GeV ILC [45]. The 
scenario with sizable singlet/doublet mixing between the and /ig states can be probed 
in the Higgs-strahlung process, e~^e~ — > Z+Higgs, where the separate Higgs states can be 
disentangled even if they are nearly degenerate in mass, as the resolution on the Higgs 
masses in this process is smaller than 100 MeV [45]. The scenario of eq. fl4.29p . in which 
there is a light CP-even Higgs particle, allows the — > h^hi decay to occur, which can 
also be probed in the Higgs-strahlung process e~^e~ h^^Z fi~^fi~bbbb, and both and 
h?, masses could be accurately determined. The singlet-like state a? could be also accessible 
in the pair production process e~^e~ /i^a? bbbb unless the mass is too large or the 
coupling A prohibitively tiny. 

4.4 Direct and indirect detection of dark matter 

For completeness, let us comment on the prospects for detecting the singlino-like LSP Xi 
(with non-singlet components of 0{X)) in astroparticle experiments. For the direct LSP 
detection via Xi-nuclei interactions, the prospects are quite dim: due to the very small values 
of A that are required from compatibility of the Higgs sector with LEP bounds, one finds 
extremely small WIMP-nucleon cross section^. We computed these cross sections using 
the recent version of the dark matter code MicrOMEGAS [47], adapted to the NMSSM. 
The result is that the relevant values for the spin-independent and spin-dependent proton 
and neutron interaction cross-sections are always below ~ 10~^^ pb, implying an expected 
number of events smaller than 0(lO~^"^~^)/day/kg for both ^^Ge and ^^^Xe nuclei, and less 
than O{10-^^^^^^)/daY/kg for ^He. 

Furthermore, since Xi has extremely small couplings to SM fermions and gauge bosons, 
the annihilation cross sections of Xi into SM particles are always extremely small. The only 
significant channel could be the annihilation process XiXi ~^ '^1 with the CP-odd singlet- 

''^In different studies [18,46], large direct detection cross sections have been obtained for the case of the 
non-constrained or semi-constrained NMSSM. The large cross sections were associated to lighter Xi and 
larger values of A, which arc not possible in the present case. 
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like state being sufficiently light to be produced on-shell. However, the a^XiXi coupling 
is proportional to k which is tiny in our case, k < 10~^. 

Thus all cross sections involving the singlino LSP are extremely small; in other words, 
the fully constrained cNMSSM can be excluded by the direct or indirect detection of a 
WIMP-like dark matter candidate. 

5 Summary and conclusions 

Due to its simplicity, namely a scale invariant superpotential and universal SUSY-breaking 
terms at the GUT scale, the cNMSSM is a very attractive supersymmetric extension of 
the standard model. We have found that for small values of the NMSSM specific Yukawa 
coupling, A < 10~^, the cNMSSM can satisfy all present phenomenological constraints 
including LEP constraints on the Higgs sector. The correct dark matter relic density is 
obtained for small mo and Aq as compared to Mi/2- Moreover, in the region M1/2 ~ 
550 — 600 GeV, the deviation of {g — 2)^ from its SM model value and the two excesses of 
Higgs-like events observed at LEP can be simultaneously explained. 

The cNMSSM sparticle spectrum allows to discriminate it from most versions of the 
MSSM: all squarks are lighter than the gluino, and due to the weakly coupled singlino-like 
LSP and the stau NLSP, the latter will show up in practically all sparticle decay cascades. 
In some regions of the parameter space, the stau lifetime can be sufficiently large, possibly 
leading to visibly displaced vertices. 

The mass of the SM-like CP-even Higgs boson is constrained to lie in the 115 - 120 GeV 
range, and for certain regions in parameter space it could strongly mix with the singlet-like 
Higgs state. In a small - different - region of the parameter space, decays of the SM-like 
CP-even Higgs boson into two singlet-like Higgs states are possibly detectable at an ILC. 

Note that the model is very predictive: as it becomes clear from Figs. [6] and [TJ the mea- 
surement of one sparticle mass (or mass difference) would allow to predict quite accurately 
the complete remaining sparticle spectrum. On the other hand, due to the singlino-like 
LSP together with small A, dark matter detection signals can rule out the present model. 

Hopefully, the cNMSSM can be tested in the near future at the LHC. To this end the 
sensitivities of the ATLAS and CMS detectors to sparticle decay cascades involving the 
stau NLSP should be thoroughly studied. 
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